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ABSTRACT 

We evaluate the effect of electrostatic screening by ions and electrons on 
low-Z thermonuclear reactions in the sun. We use a mean field formalism and 
calculate the electron density of the screening cloud using the appropriate 
density matrix equation of quantum statistical mechanics. Because of well 
understood physical effects that are included for the first time in our treatment, 
the calculated enhancement of reaction rates does not agree with the frequently 
used interpolation formulae. Our result does agree, within small uncertainties, 
with Salpeter's weak screening formula. If weak screening is used instead of the 
commonly employed screening prescription of Graboske et al., the predicted 8 B 
neutrino flux is increased by 7% and the predicted chlorine rate is increased by 
0.4 SNU. 



Subject headings: nuclear reactions 



1. Introduction 

In recent years, an increasing amount of attention has been devoted to calculating 
more accurately the effects on the rates of solar fusion reactions of electrostatic screening in 
the solar plasma (Carraro, Schafer, & Koonin 1988; Johnston, Kolbe, Koonin, & Langanke 
1992; Bahcall & Pinsonneault 1992; Shoppa, Koonin, Langanke, & Seki 1993; Dzitko, 
Turck-Chieze, Delbourgo-Salvador, & Lagrange 1995; Ricci, Degl'Innocenti, & Fiorentini 
1995; Gruzinov & Bahcall 1997; Brown k Sawyer 1997; Briiggen & Gough 1997). All 
of these discussions take as their starting point the classical analysis by Salpeter (1954). 
The primary reason for making more precise calculations is that nuclear fusion reactions 
produce the solar neutrino fluxes (see |http:/ /www.sns.ias.edu^jnb| ). The neutrino fluxes 



are being observed with a number of large new detectors that are expected to yield flux 
measurements of high accuracy (of the order of a few percent or better, see Bahcall et al. 
1995, and, for more details, Totsuka 1996; McDonald 1994; Arpesella et al. 1992). 
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In this paper, we calculate for the first time the electron density in the vicinity of the 
fusing nuclei using the partial differential equation for the density matrix that is derived 
in quantum statistical mechanics. In previous treatments of screening that attempted 
to go beyond the linear regime, the electron density near the nucleus was either taken 
to be-without quantitative justification-the unperturbed value, n e (oo) (Mitler 1977; 
Dzitko, Turck-Chieze, Delbourgo-Salvador, & Lagrange 1995) or left as a free parameter 
(Ricci et al. 1995) or the electrons were assumed to be completely degenerate (Graboske, 
DeWitt, Grossman, & Cooper 1973). We calculate screening corrections in a mean field 
approximation; we numerically solve the nonlinear Poisson-Boltzmann equation for a 
mixture of electrons and ions. The electron density distribution calculated from the density 
matrix equation is included self-consistently and iteratively in the mean field equation. 

Our results represent both an improvement on, and a simplification of, the description 
of nuclear fusion used in many solar evolution codes. 

For simple physical reasons, our results differ from the interpolation formulae that are 
currently used to describe reaction rates in the Sun (Salpeter & Van Horn 1969; Graboske 
et al. 1973), and the numerical calculations of Dzitko, Turck-Chieze, Delbourgo-Salvador, 
& Lagrange (1995). 

Interpolation formulae describe a transition between Salpeter's weak screening, which 
is due to both electrons and ions, and strong screening, for which only ions are effective. At 
the high densities relevant for strong screening, electrons are fully degenerate. The solar 
core, however, is only weakly degenerate, and the effects of degeneracy are already included 
into the Debye radius, Rd (which is increased 2% by electron degeneracy, cf. Eq. ([|) of 
the present paper or Eq. (25) of Salpeter 1954). Therefore, the interpolation formulae in 
use underestimate the electron contribution to screening and give reaction rates lower than 
ours. 

The numerical procedures of Dzitko, Turck-Chieze, Delbourgo-Salvador, & Lagrange 
(1995) and Mitler (1977) predict reaction rates that are too slow for heavy ions because 
they assumed that the electron charge density near a screened nucleus is the unperturbed 
value, en e (oo). This assumption seriously underestimates the charge density near heavy 
ions. For example, it is known that a screened beryllium nucleus under solar interior 
conditions has charge density near the nucleus ~ — 3.85en e (oo) (Gruzinov & Bahcall 1997, 
Brown & Sawyer 1997; all quantum mechanical calculations give similar results, see Bahcall 
1962 and Iben, Kalata, & Schwartz 1967). 

This paper is organized as follows. In §2 we review the basic concepts and in §3 
we relate the electrostatic energy to the screening enhancement using the free energy. 
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We describe the calculations in §4 and summarize the numerical results in §5. In §6 we 
summarize our main results and present the conclusions regarding solar neutrino fluxes. 
The Appendix evaluates a quantum correction to the kinetic energy of thermal electrons in 
the electrostatic field of a screened nucleus. 



2. Enhancement of Fusion Rates 

The solar core plasma is dense enough that it noticeably enhances fusion rates as 
compared to the rates in a rarefied plasma of the same temperature. As explained by 
Salpeter (1954), the rate of a fusion of two nuclei of charges Z\ and Z 2 is increased by a 
factor 

/ = expA, (1) 



where 



Here Rd is the Debye radius, 



e 2 

A = Z X Z 2 . (2) 
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Here (3 = 1/T, n is the baryon density, X iy Z^ and Ai are, respectively, the mass 
fraction, the nuclear charge, and the atomic weight of ions of type i. The quantity 
/'// ~ 0.92 accounts for electron degeneracy. Equation (f|) is the same as Eq. (25) of 
Salpeter (1954). In what follows, we will make use of a simplified expression for £, 

Cs im ple = {(l-r/2)W+l} 1/2 , (5) 

in which the plasma is assumed to consist only of hydrogen and helium (Y is the helium 
abundance by mass). The approximation of considering only a hydrogen and helium plasma 
rather than the full solar composition (cf. Grevesse & Noels 1993) causes an error of less 
than 0.5% in computing solar fusion rates. This error is completely unimportant for our 
purpose of estimating the ratio of the total screening to the weak screening value given by 
Eqs. (D-®. 

The enhancement of fusion rates due to screening depends only very weakly upon 
location in the solar interior (cf. Ricci et al. 1995), because the primary effect of screening 
is proportional to p/T 3 (cf. Eq. |1]-^), which is approximately constant in the solar interior 
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(cf. Eq. 41 of Bahcall et al. 1982). The plasma parameters at a characteristic radius in the 
solar interior, R/R & = 0.06, are (Bahcall & Pinsonneault 1995) Rr, = 0.46 and T = 47 in 
atomic units (m e = h = e = 1). In units that are more common in astronomical discussions, 
the temperature, Tq, in millions of degrees is Tq = 0.32T = 15 and the Debye radius in cm 
is R D = 0.46 x 5.3 x 10~ 9 = 2 x 10" 9 cm. 

Consider an important example: Z\Z 2 = 4 for the solar fusion reactions 3 He( 4 He, 7) 7 Be 
and 7 Be(p,7) 8 B. Eq. (2) yields A = 0.19. According to Eq. (1), the calculated rates of 
these fusion reactions are then 21% faster than they would be if screening were neglected. 

Eq. (1) is only valid to first order in A. Nonlinearities in the electrostatic screening 
interactions might naively be expected to produce corrections ~ A 2 , i.e., of order 4% for a 
Z\Z 2 =4 reaction. In the following sections, we calculate corrections to the Salpeter weak 
screening formula, Eq. (|l|), and find that the numerical corrections are always significantly 
smaller than A 2 . 



3. Enhancement Factors and Free Energy 

Salpeter 's formula (Eqs. (1),(2)) can be derived as follows. The screened potential near 
the nucleus Z\ in the Debye-Hiickel approximation is 

ps — - — . (6) 

r r Ro 

The potential shift Z^jRu increases the probability that the charge Z 2 comes close to Z\ 
by the Boltzmann factor e A , A = ^Z\Z 2 j Rd- 

Unfortunately, this clear derivation can not be used if we go beyond the Debye-Hiickel 
approximation and include nonlinear screening effects. Given a numerically calculated 
potential around the charge Z\, 0i(r) , we can not assume that the enhancement factor is 
equal to e A with A = (3Z 2 x (Z\/r — <f>i(r))\ r= o. This is already obvious from the asymmetry 
of this expression under the 1-2 permutation; <pi is not just proportional to Z\ for nonlinear 
screening. 

In the more general case considered here, the enhancement of fusion rates due to 
screening can be calculated in terms of an expression involving the free energy of a screened 
charge Z, F(Z). In terms of free energy, the enhancement factor is simply (DeWitt, 
Graboske, & Cooper 1973) e A with 



-I3F(Z 1 + Z 2 ) + (3F(Z l ) + (3F(Z 2 ), 



(7) 
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which is a manifestly symmetrical expression. Equation(^) expresses the thermodynamic 
relation that at constant temperature (which is relevant when considering solar fusion 
reactions) 5F = — SW, where 5W is the work done by the plasma on the fusing ions. 
The extra work performed by the plasma due to screening is positive, pushing the fusing 
ions closer together. For a given relative kinetic energy when the ions fuse, the initial 
kinetic energy is lower by 5W than in the absence of screening. Therefore, the probability 
of the fusing configuration is increased, i.e., the reaction rates are faster, by a factor 
exp(j36W/T) = exp(-(36F) = exp(A). 

The free energy can be calculated in terms of electrostatic energy using the 
thermodynamic formula 

(3F= [ d/3'U. (8) 
Jo 

The lower limit in the integral in Eq. @ is chosen so that at high-temperature (small 
(3) F goes to zero as P 1 ^ 2 , as implied by Debye theory (see discussion below). The total 
electrostatic energy including the self-energy is 

1 



U tot = ~J d 6 r(j)(r)p(r). (9) 

The self-energy of the charges cancels out in performing the difference indicated by 
Eq. (0). The fusing nuclei are well separated whenever screening is relevant; their combined 
self-energies are the same in the fusing state as the sum of the self-energies in the initial 
(infinitely separated) state. Most of the acceleration of the fusing nuclei occurs at distances 
larger than O.IRd, which is four orders of magnitude larger than nuclear radii. Therefore, 
the relevant self-energy for the calculation of enhancement factors due to screening does not 
include the self-energies and is 

U = ^5<P(0) + lfd 3 r<P(r)5p(r), (10) 

where 5<fi = <fi — Z/r and Sp = p — Z8(r). 

In the Debye-Hiickel approximation these expressions reproduce Salpeter's formula 
(Briiggen & Gough 1997). In the Debye-Hiickel approximation, = (Z/r)e~ r ' RD , 
4ir5p = -<j)/R 2 D , 5(f)(0) = -Z/R D , and Eq. (0) gives U = -(3/4)Z 2 /R D . Since 
R D ~ /T 1 / 2 , Eq. © gives (3F = -(l/2)f3Z 2 / R D . Then Eq. (|) gives Eq. (|). 



- 6- 



4. Calculations 

In the mean field approximation, electrostatic screening of a charge Z is described by 
the Poisson-Boltzmann equation 

V 2 = 4vm{(l - -)e w - (1 - Yy-M - -e~ 2 ^}, (11) 

where the terms on the right hand side represent, respectively, screening by electrons, 
protons and alphas. The boundary condition is —>■ Z/r for r — > 0. In the non-linear 
regime, one cannot solve the Poisson-Boltzmann equation as written. Classical electrons 
recombine, which corresponds formally to the divergence of the classical Boltzmann 
factor near the nucleus. This problem does not arise in previous solutions of the 
Poisson-Boltzmann equation which were carried out in the linear regime corresponding to 
weak screening. 

Quantum statistical mechanics must be used for calculating terms beyond the weak 
screening approximation. Fortunately, electron degeneracy makes only a small correction, 
less than 2% in the Debye radius, see Eq. |], and therefore less than 1% for all cases in the 
reaction rates. Hence, a distinguishable particles approximation can be employed^ We use 
a numerical code that solves the density matrix equation for the density of electrons near 
the nucleus. The code, which was developed following the discussion of Feynman (1990), is 
described in Gruzinov & Bahcall (1997). 

The average electron density can be calculated by solving the density matrix equation 
(e.g. Feynman 1990) 

V={^V 2 + 0(r)}p, (12) 

with the initial condition 

p(r,/3 = 0) = 5 (3) (r). (13) 

Since Eq. ( |T2"D appropriately describes the quantum statistical mechanical effects, the 
solution for the density matrix converges everywhere despite the divergence of the classical 
potential at r = 0. Another great advantages of the density matrix formulation is that the 
character of the states in the plasma does not have to be specified and therefore difficult 
questions concerning the existence or non-existence of bound states are finessed. The 



1 Brown and Sawyer (1997) calculated electron densities using both Fermi-Dirac and Maxwell-Boltzmann 
statistics. Degeneracy effects on the value of the central electron density were of order 10% for Z = 6. 
We shall show in the course of this paper that changing the central electron density by almost an order of 
magnitude does not significantly change the rate of nuclear fusion reactions. Therefore, the small fractional 
change in the central electron density due to using different statistics is not important for our purposes. 
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enhancement of the electron density to be used in the Poisson-Boltzmann equation instead 
of the Boltzmann factor is the solution of Eq. ([12]) for the nuclear charge of Z divided 
by the solution for Z = 0. The solution for the Z = case can be obtained analytically and 

is po(/3) = (27T/3)" 3 / 2 . 

As described in Gruzinov & Bahcall (1997), the diffusion with multiplication problem, 
Eq. fljjP, can be solved easily by direct three-dimensional numerical simulations for solar 
conditions, because the inverse temperature (3 is small (~ 0.02), and the diffusive trajectory 
stays close to the origin. The mesh size and the regularization procedure were the same as 
in our previous work. 

Numerically, we start with an initial guess that <f>(r) = (Z/r)e~ r ^ RD everywhere and 
then calculate the electron density using Eq. (|i"2"[ ) for all r < 0.4. The particular value of 
r = 0.4 (~ Rd) is not important. For all r ^ 0.2, our density matrix code simply reproduces 
the Boltzmann distribution factor n(r) = n(oo)e^. We use the calculated electron density 
at r < 0.4 to solve Eq. (|Tl|) numerically for all r. We then obtain a new potential </>(r). We 
use this potential to calculate the electron density at r < 0.4 using Eq. QT2j) and repeat the 
procedure. The procedure converges quickly, after one to three iterations. 

The electrostatic energy was calculated from Eq. (0). The calculation was repeated 
at higher temperatures for the purpose of estimating the free energy using Eq. (|j). 

Quantum statistical mechanics implies the existence of an effect that we believe has 
not been previously considered in the context of fusion reaction rates. The kinetic energy 
of electrons in the electrostatic field of the nucleus is no longer |T per electron. Indeed, 
the kinetic energy of electrons is increased. In the low-temperature limit this effect is the 
familiar zero-point oscillations. In the high temperature limit the effect is more subtle, but 
it can be calculated analytically. In the Appendix, we calculate the quantum statistical 
mechanics corrections to the electron kinetic energy and the resulting correction to the free 
energy. 

5. Numerical Results 

Table 1 gives the numerical results for: (i) corrections to the Debye-Hiickel electrostatic 
energy, (ii) corrections to the free energy due to the changed electrostatic energy, (iii) 
corrections to the kinetic energy of electrons, (iv) corrections to the free energy due to the 
changed kinetic energy of electrons, and (v) the total correction to free energy. 

Figure 1 explains the sources of different corrections to the Debye-Hiickel approximation 
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of screening, (i) At large distances (small 0), the plasma response is suppressed due to helium 
ions. To see this, expand Eq. (|TT|) up to the second order in 0: V 2 = 0(1 — w 13(f)) / R 2 D , 
where w = 3Y/(8 — 2Y). (ii) At small distances, the plasma response is suppressed due to 
the fuzziness of quantum electrons, which is expressed by the density matrix Eq. flT2"|). (iii) 
At intermediate radii, the plasma response can be enhanced just because > [3(f). 

The second column of Table 2 shows the corrections, —5 A, to reaction rates calculated in 
this paper relative to Salpeter's weak screening rates. For example, the correction to the rate 
of the 7 Be(p, 7 ) 8 B reaction is 8A = -f35F(5) + f35F(4) + f35F(l) = -5.2% + 3.5% + 0.2% = 
—1.5%. This means that the reaction is only 1.5% slower in the Sun than predicted by the 
Salpeter formula. Table 2 also compares our corrections with those predicted by Graboske, 
DeWitt, Grossman, & Cooper (1973) (GDGC), by Salpeter & Van Horn (1969) (SVH), and 
by Dzitko, Turck-Chieze, Delbourgo-Salvador, & Lagrange (1995) (DTDL). 

Our corrections are typically an order of magnitude smaller than the corrections 
calculated by GDGC (cf. columns two and three of Table 2). The intermediate screening 
prescription of GDGC (their Table 4, page 465) uses the intermediate screening formula 
from DeWitt, Graboske, & Cooper (1973) (their Eq. (70), page 455); the DeWitt et al. 
formula was obtained as an illustration assuming completely degenerate electrons, which 
is inappropriate for the solar interior. The GDGC intermediate screening prescription 
underestimates the enhancement of fusion reactions by a factor, exp(<5A G DGC — <5A G b), 
which varies from about 7% for the important 3 He( 4 He, 7) 7 Be and 7 Be(p,7) 8 B reactions to 
about 16% for the 14 iV(p,7) 15 reaction (cf. Table II of Ricci et al. 1995). 

The discrepancies between our results and those of Dzitko et al. (1995) only become 
large when relatively heavy nuclei are involved. In this case, the electron density in the 
vicinity of the fusing nuclei is much larger than the value, n e (oo), assumed by Dzitko et al. 
Table II shows, for example, that when the electron density is calculated from the density 
matrix equation the value for — 5A is only 0.8% instead of the Dzitko et al. value of 6.3% 
(cf. columns two and five of Table 2). 

For heavier nuclei like nitrogen, the large classical enhancement of electron density 
near the nucleus competes with the smearing effect due to quantum fuzziness, resulting in 
a net correction that is smaller than for the lighter nuclei (see Figure 1). 

6. Summary and Conclusion 

We use the density matrix equation to determine from quantum statistical mechanics 
the electron density in the near vicinity of the fusing nuclei. Our treatment is the first to 
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describe properly the electron density in screening calculations that are appropriate for solar 
interior conditions. Previously, the lack of understanding of what to use for the electron 
density near the fusing nuclei has been the principal cause for uncertainty in estimating 
non-linear corrections to screening calculations (see, e.g., Ricci et al. 1995 and references 
therein). 

The non-linear corrections that we calculate to the Salpeter weak screening formulae, 
Eq. (^j)-Eq. (^), are, for solar conditions, ~ 1% for all the important nuclear fusion 
reactions. The principal uncertainty in our calculations is caused by thermal fluctuations, 
which are not included in the present treatment. For the analogous case of electron capture, 
thermal fluctuations affect the average rate by < 1% (Gruzinov and Bahcall 1997). Since 
the non-linear effects calculated in the present paper are small and of the same order as the 
effects of fluctuations that occur in the electron capture problem, we recommend using the 
Salpeter weak screening formula for solar fusion rates. 

What difference do the present results make for the solar neutrino problem? This 
question is answered by Table 3 of Bahcall & Pinsonneault (1992). Keeping all other input 
data constant, the weak screening approximation gives, relative to the Graboske et al. 
prescription, a 0.4 SNU larger result in the chlorine (Homestake) experiment, a 2 SNU 
increase in the gallium experiments, and a 7% larger 8 B neutrino flux (measured in the 
Super Kamiokande, Kamiokande, and SNO experiments). The Graboske et al. prescription 
was used previously by Bahcall and Pinsonneault and in many other stellar evolution codes 
(cf. Ricci et al. 1995). 

An error in the screening enhancement is equivalent to an error in the low energy cross 
section factor. Therefore, one can use the well known power law dependences of the neutrino 
fluxes on cross section factors (Bahcall 1989) to estimate the uncertainties introduced by 
inaccuracies in the screening calculations. A 1% uncertainty in the screening calculation 
causes an ~ 1% uncertainty in the predicted 8 B neutrino flux and a smaller uncertainty 
for other fluxes in the pp chain. For the crucial 8 B neutrino flux, the uncertainty in the 
measurement of the low energy cross section factor for the 7 -Be(p,7) 8 5 reaction causes a 
much larger uncertainty, > 10% (see Bahcall and Pinsonneault 1995). 

The non-linear effects in ion and electron screening that are evaluated in this paper 
cause differences in the solar model neutrino fluxes that are small compared to the 
order-of-unity differences between the rates measured in solar neutrino experiments and the 
fluxes predicted by standard models (assuming nothing happens to the neutrinos after they 
are created). 



This work was supported by NSF PHY-9513835. We are grateful to S. Turck-Chieze for 
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valuable discussions that first directed our attention to the problem of the large apparent 
corrections implied by the Graboske et al. prescription and stimulating comments on a 
draft of this manuscript. We are grateful to M. Briiggen for a valuable discussion. 

A. Quantum Corrections to the 3/2T Kinetic Energy Per Particle Rule 

In classical statistical mechanics, the kinetic energy of particles, interacting or 
noninteracting, in an external potential or in free space, is §T per particle. In quantum 
statistical mechanics, the kinetic energy at a given temperature depends on the external 
potential. This is obvious in the low-temperature limit: the kinetic energy of the ground 
state is positive if the external potential is not zero (the zero-point oscillations). 

Thermal electrons in an electrostatic field of a nucleus have kinetic energy larger than 
|T. The effect depends on Z, and reduces the reaction rates (as compared to Salpeter's 
weak screening rates). The correction to kinetic energy can be calculated if the diagonal of 
the density matrix (e.g., Feynman 1990) p(r, f3) = p(r,r,(3) is known, 

SK = n e (27r(3 f 2 J d 3 r{-dpp - (^/T 1 + V)p}, (Al) 

that is the correction to the kinetic energy is the total energy minus the unperturbed kinetic 
energy, |T, minus the potential energy, V. In classical statistical mechanics 

p = (2 7 r/?)- 3 / 2 e- /3y , (A2) 

and Eq. (Al) gives 5K = 0. 

We calculated the density matrix of electrons analytically and used Eq. (Al) to 
calculate the kinetic energy correction for V = — fe~ r / RD assuming (3 <C 1, Z ~ few, 
Rd > These conditions are satisfied in the solar interior where (3 ss 0.02 and Rd ~ 0.5. 
Two different approaches were used at distances from the nucleus greater than the de 
Broglie wavelength, /3 1 / 2 , and at distances smaller than the Debye radius R D . These two 
approaches are explained below. 

A.l. r ^> [3 1 / 2 : High- Temperature Expansion 

Thermal electrons have "a characteristic size" ~ (3 1 / 2 . If the potential energy does not 
change by much over this distance (which in our case is true for r > /9 1//2 ), the density 
matrix is approximately given by Eq. (A2) with small corrections. The corrections are due 
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to the fact that a fuzzy thermal electron samples potential not only at a given point but in 
the /3 1//2 -vicinity of the given point. 

Let (x,y,z) be a small deviation of coordinates from (r, 0,0). Potential energy is, up 
to the second order, 

W = V'x + \v»f + \^f+ 1 ^ (A3) 

where prime denotes the r-derivative, and we assume that V is spherically symmetrical. The 
path integral giving the density matrix (e.g., Feynman 1990, chapter 3) is Gaussian and can 
be calculated. In fact, the answer can be constructed without the actual calculation from 
the known density matrix of the linear harmonic oscillator (e.g., Feynman 1990, chapter 2). 
It reads 

p = (27r/3)- 3/2 e-^{l + ^V' 2 - \-B\V" + V)}. (A4) 

24 12 r 

The kinetic energy correction is given by Eq. (Al) 

5K = n e J 47ir 2 dre^ v {-^ 2 V' 2 + ^(3(V" + *V')}. (A5) 

In our case Eq.(A4) is valid only at r ^ P 1 ^ 2 , but if potential energy V were smooth at all 
r, we could have integrated the last term by parts 

SK = ^n e [3 2 J An^dre-^V' 2 , (A6) 

showing that kinetic energy correction is positive in the high-temperature limit. In our 
calculation we used Eq. (A4) at r > r , and results from the next section were used at 
r < ro . The final answer does not depend on the choice of r as long as Rd ^ r ^ /3 1 / 2 . 



A. 2. r <C R D : Hydrogenic Density Matrix 

At distances from the screened nucleus r <C Rd, the potential energy is 

1 =-7 exp( -^ )K -7 + iV (A7) 

The only effect of the constant correction is to lower electron density by the Boltzmann 
factor e -$ z l R D _ The density matrix in the Coulomb potential can be obtained from 
hydrogenic eigenstates. 

The kinetic energy correction is 

5K = n e e~^ {2n(3f /2 J d 3 r{-dp P - (-/T 1 + V)p}. (A8) 
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The density matrix is 

00 9/ -I- 1 00 a r°° dk h* 2 

p(r,P)=E-^r{E\Rn l (r)fe^+ 5^)1 V^}. (A9) 

«=0 ^ n=l J ° Z7F 

Here the bound states of hydrogen are (e.g. Landau & Lifshitz 1977) 

Rm(r) = nl+2{ l - 1} , { } 1/2 ^' r/n n-n + 1 + 1,21 + 2, 2r/n), (A10) 

where F is the confluent hypergeometric function. The continuum states are 

R ki {r) = 2ke 7r/2k \T(l + 1 - % -)\{2kr) 1 e~ ikr F Ct + 1 + 1,21 + 2, 2ikr), (All) 

K K 

and for Z 7^ 1 we scale r — > Zr, j3 — > Z 2 /3. 

We used these formulae to calculate the kinetic energy shift at small r. Results of 
this subsection match the high-temperature results if /3 1 / 2 < r < R^. We repeated the 
calculation at smaller (3 to obtain the free energy shift due to the quantum correction to 
kinetic energy of electrons, 

(35 F= [ dp'SK(p'). (A12) 
Jo 

Results are shown in Table I. 



REFERENCES 

Arpesella, C, et al. 1992, BOREXINO Proposal, Vols. 1 and 2, ed. G. Bellini, et al. 
(Milano: Univ. of Milano). 

Bahcall, J. N. 1962, Phys. 128, 1297. 

Bahcall, J. N., Lande, K., Lanou, R. E., Jr., Learned, J., Robertson, R. G. H., & 
Wolfenstein, L. 1995, Nature, 375, 29 

Bahcall, J. N., Huebner, W. F., Lubow, S. H., Parker, P. D., and Ulrich, R. K., 1982, Rev. 
Mod. Phys. 54, 767. 

Bahcall, J. N. 1989, Neutrino Astrophysics (Cambridge: Cambridge University Press) 
Bahcall, J. N., & Pinsonneault, M. H. 1992, Rev. Mod. Phys., 64, 885 
Bahcall, J. N., & Pinsonneault, M. H. 1995, Rev. Mod. Phys., 67, 781 
Brown, L. S., & Sawyer, R. F. 1997, ApJ, 489, 968 
Brviggen, M., & Gough, D. O. 1997, ApJ, 488, 867 



-13- 



DeWitt, H. E., Graboske, H. C, & Cooper, M. S. 1973, ApJ, 181, 439 

Dzitko, H., Turck-Chieze, S., Delbourgo-Salvador, P., & Lagrange, C. 1995, ApJ, 447, 428 

Carraro, C, Schafer, A., & Koonin, S. E. 1988, ApJ, 331, 565 

Feynman, R. P. 1990, Statistical Mechanics (Reading: Addiso n- Wesley) , Chapters 2 and 3 

Graboske, H. C, DeWitt, H. E., Grossman, A. S., & Cooper, M. S. 1973, ApJ, 181, 457 

Grevesse, N., & Noels, A. 1993, Origin and Evolution of the Elements, ed. N. Prantzos, E. 
Vangioni-Flam, & M. Casse (Cambridge: Cambridge Univ. Press), 15 

Gruzinov, A. V., & Bahcall, J. N. 1997, ApJ, 490, 437 

Iben, L, Jr., Kalata, K., & Schwartz, J. 1967, ApJ, 150, 1001. 

Johnston, C. W., Kolbe, E., Koonin, S. E., & Langanke, K. 1992, ApJ, 392, 320 

Landau, L. D., & Lifshitz, E. M. 1977, Quantum Mechanics, Third Edition (Oxford: 
Pergamon Press) 

McDonald, A. B. 1994, Proceedings of the 9th Lake Louise Winter Institute, ed. A. 
Astbury, et al. (Singapore: World Scientific), 1 

Mitler, H. E. 1977, ApJ, 212, 513 

Ricci, B., Degl'Innocenti, S., & Fiorentini, G. 1995, Phys. Rev. C, 52, 1095 
Salpeter, E. E. 1954, Australian J. Phys., 7, 373 
Salpeter, E. E., & Van Horn, H. M. 1969, ApJ, 155, 183 

Shoppa, T. D., Koonin, S. E., Langanke, K., & Seki, R. 1993, Phys. Rev. C, 48, 837 
Totsuka, Y. 1996, Nucl. Phys. B (Proc. Suppl.), 48, 547 



This preprint was prepared with the AAS IATeX macros v4.0. 



-14- 



Table 1. Electrostatic, Kinetic and Free Energy Corrections 



z 


1 


2 


4 


5 


7 


8 


06U, % 


0.34 


1.6 


6.4 


9.2 


11.2 


7.6 


(3SFu,% 


0.1 


0.6 


2.7 


3.9 


5.7 


5.2 


05K, % 


0.22 


0.57 


1.9 


3.2 


8.1 


12.6 


0SF K ,% 


0.1 


0.3 


0.8 


1.3 


2.9 


4.4 


0SF, % 


0.2 


0.9 


3.5 


5.2 


8.6 


9.6 



Note. — The symbols represent: nuclear charge Z, corrections to: (i) the electrostatic energy normalized 
to temperature 0SU, (ii) the free energy due to increased electrostatic energy 0SFu 7 (iii) the kinetic energy 
of electrons 06K, (iv) the free energy due to increased kinetic energy 0SFk, (v) the total free energy 08F. 
The plasma parameters are taken from the solar model of Bahcall & Pinsonneault 1995 at the representative 
point R/Rq = 0.06. 



Table 2. Reaction Rate Corrections 



Reaction 


GB 


GDGC 


SVH 


DTDL 


P+P 


0.5 


0.0 


0.5 


0.2 


3 He+ 4 He 


1.7 


8.2 


2.4 


1.8 


p+ 7 Be 


1.5 


8.5 


2.6 


2.3 


p+ 14 N 


0.8 


15.2 


6.3 


6.3 



Note. — Corrections to weakly screened reaction rates. Nuclear fusion reactions in the sun are enhanced 
by a factor exp(A + SA), where A is given by Salpeter's expression, which is Eq. |[ The table shows the 
corrections, — 5A in percent, calculated in this paper (GB), by Graboske, DeWitt, Grossman, & Cooper 
(1973) (GDGC), by Salpeter & Van Horn (1969) (SVH), and by Dzitko, Turck-Chieze, Delbourgo-Salvador, 
& Lagrange (1995) (DTDL). The corrections refer to a representative point R/Rq = 0.06. 
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Fig. 1. — Screening of two test charges. The induced charge density, p, normalized to the 
Debye-Huckel charge density, pdh = 0/(4vri?|)), is shown as a function the electrostatic 
potential, 0. For small (large distances from the screened nucleus), p is given by the 
classical Boltzmann formula (see the right hand side of Eq. |TTJ ). In this region, p is smaller 
than pdh due to the presence of helium ions. At large (p (close to the screened nucleus), 
p is much smaller than pdh due to the quantum fuzziness embodied in the density matrix 
equation (Eq. [T2|). At intermediate distances, p can be large than p DH (because exp0 > (f> 
). For Z = 8 the plasma response is larger than linear at intermediate distances. For Z = 5 
the plasma response is always smaller than linear. 



